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ABSTRACT 

With a goal toward deriving the physical conditions in external galaxies, we present a survey of 
the formaldehyde emission in a sample of starburst systems. By extending a technique used to derive 
the spatial density in star formation regions in our own Galaxy, we show how the relative intensity of 
the lio — In and 2n — 2i 2 K-doublet transitions of H2CO can provide an accurate densitometer for 
the active star formation environments found in starburst galaxies. Relying upon an assumed kinetic 
temperature and co-spatial emission and absorption from both H 2 CO transitions, our technique is 
applied to a sample of nineteen IR-bright galaxies which exhibit various forms of starburst activity. 
In the five galaxies of our sample where both H2CO transitions were detected we have derived spatial 
densities. We also use H 2 CO to estimate the dense gas mass in our starburst galaxy sample, finding 
similar mass estimates for the dense gas forming stars in these objects as derived using other dense 
gas tracers. A related trend can be seen when one compares Lm to our derived n(.ff 2 ) for the five 
galaxies within which we have derived spatial densities. Even though our number statistics are small, 
there appears to be a trend toward higher spatial density for galaxies with higher infrared luminosity. 
This is likely another representation of the Lm-Mdense correlation. 
Subject headings: galaxies: starbursts, ISM: molecules 



1. INTRODUCTION 
Studies of the distribution of Carb on Monoxide (CO) 



emissi on in external galaxies (cf. Young & Scoville 



(1991)) have pointed to the presence of large quantities 
of molecular material in these systems. These studies 
have yielded a detailed picture of the molecular mass 
in many external galaxies. But, because emission from 
the abundant CO molecule is generally dominated by ra- 
diative transfer effects, such as high optical depth, it is 
not a reliable monitor of the physical conditions, such as 
spatial density and kinetic temperature, quantities neces- 
sary to assess the possibility of star formation. Emission 
from less-abundant, higher-dipole moment molecules are 
better-suited to the task of deriving the spatial density 
and kinetic temperature of the dense gas in our and exter- 
nal galaxies. For this reason, emission line studies from 
a variety of molec ules have been made toward m ainly 
nearby galaxies (s ee Mauersberger & Henkel ( 1989 ) (CS), 
Gao fc Solomon] j2fMa| (HCJN ), |Nguyen- G-kleu i ^L«i 
992| (HCO^ [Mj^er_ sberger et oZ] fll99o| and"^ 



Lurncr 



(NH 3 





auersberger et al. 



auersberger "| 1991 ) for a 
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review) . 

The most extensive sets of measurements of molecular 
line emission in external gal axies has been don e using the 
J=l-0 transitions of CO dHelfer et aZ.||2003| and HCN 
dGao & Solomonl 12004a I. 



20031 

=1-0 transitions 
of CO and HC1N are good tracers of the more generally 
distributed and the denser gas, respectively, but do not 
provide comprehensive information about the individual 
physical conditions of the dense, potentially star-forming 
gas, another molecule must be observed for this purpose. 
Formaldehyde (H2CO) has proven to be a reliable den- 
sity and kinetic temperature probe in Galactic molecular 
clouds. Existing measurements of the H 2 CO lip— I n and 
2n — 2i2 emission in a w i de va ri ety of galax i es by Baan 
et a?.|(11986| |Baan et gj] ( |1990[ ), [Baan et al.\ ( |1993[ ), and 



Araya et ai \ ( 2004[ ) have mamly concentrated on mea- 
surements of the lio — In transition. One of our goals 
with the present study was to obtain a uniform set of 
measurements of both K-doublet transitions with which 
the physical conditions, specifically the spatial density, 
in the extragalactic context could be derived. Using the 
unique density selectivity of the K-doublet transitions of 
H 2 CO we have measured the spatial density in a sample 
of galaxies exhibiting starburst phenomena and/or high 
infrared luminosity. In Sj2]we discuss the specific proper- 
ties of the H 2 CO molecule which make it a good probe of 
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Fig. 1. — H2CO K-doublet excitation temperatures as a function 
of molecular hydrogen density at an ortho-H2CO column density 
per velocity gradient of 10 1 cm - 2 /km s" 1 and kinetic temper- 
ature of 40 K. Note how the J<5 K-doublet transitions go into 
absorption for n(H2) < 10 5,5 cm" 3 . 



spatial density. <|3] presents our observation summary; £|4] 
our H2CO, OH, Hllla, and continuum emission mea- 
surement results; Sj5] analyses of our H 2 CO, OH, and 
Hllla measurements, including Large Velocity Gradi- 
ent (LVG) model fits to and dense gas mass calculations 
based on our H2CO measurements. 

2. FORMALDEHYDE AS A HIGH DENSITY PROBE 

Formaldehyde (H 2 CO) is a proven tracer of the high 
density environs of molecular clouds. It is ubiquitous: 
H 2 CO is associated with 80% of the HII regions surveyed 
by Downes et al. (19801, and possesses a large number 
of observationally accessible transitions from centimeter 
to far-infrared wavelengths. Because H 2 CO is a slightly 
asymmetric rotor molecule, most rotational energy lev- 
els are split by this asymmetry into two energy levels. 
Therefore, the energy levels must be designated by a to- 
tal angular momentum quantum number, J, the projec- 
tion of J along the symmetry axis for a limiting pro- 
late symmetric top, K_i, and the projection of J along 
the symmetry axis for a limiting oblate symmetric top, 
K + i. This splitting leads to two basic types of transi- 
tions: the high-frequency AJ = 1, AK_ X = 0, AK +1 = 
— 1 "P-branch" transitions and the lower- frequency A J 
= 0, AK_i = 0, AK+i = ±1 "Q-branch" transitions, 
popularly known as the "K-doub l et" tr ansitions (see dis- 
cussion in Mangum & Wootten ( 1993 1). The P-branch 
transitions are only seen in emission in regions where 
n(H 2 ) > 10 4 cm" 3 . The excitation of the K-doublet 
transitions, though, is not so simple. For n(H 2 ) < 10 5 ' 5 
cm~ 3 , the lower energy states of the l i0 — In through 
5i4 — 5i5 K-doublet transitions be come overpopulated 



due to a collisional se lection effect (Evans et al. (19751 
Garrison et aD (19751). This overpopulation cools the 
J < 5 K-doublets to excitation temperatures lower than 
that of the cosmic microwave background, causing them 
to appear in absorption. For n(H 2 ) > 10 5 5 cm -3 , this 
collisional pump is quenched and the J < 5 K-doublets 
are then seen in emission over a wide range of kinetic 
temperatures and abundances (see Figure n§. 



3. OBSERVATIONS 

The measurements reported here were made using the 
National Radio Astronomy Observatory (NRAO 1 ) Green 
Bank Telescope (GBT) during the period 2006/06/07- 
14. Single-pointing measurements were obtained of the 
lio - In (4.829660 GHz) and 2 n -2 12 (14.488479 GHz) 
K-doublet transitions of H 2 CO, the Hllla radio re- 
combination line (RRL) at 4.744183 GHz, and two of 
the three (F = 1—0 and 1 — 1) hyperfine satellite 
lines of the rotationally excited 2 n!/ 2 J =1/2 state of 
OH at 4.750 656(3) and 4.765562(3) GHz ( |Radfordp68 



Lovas 1992[) toward a sample of 19 infrared-luminous 
and/or star burst galaxies (Table [T]). Our source sam- 
ple was chosen to represent bo th galaxies with m e asured 
H 2 CO emis sion or absorption (Baan et al. (19931, Araya 
et al. (2004)) and other galaxies with substantial m olec- 
ular emission dedu ced from HCN measurements (Gao 
fc Solomonl ( |2004a| )). The single-beam H 2 CO li - T~n7 
Hllla, and OH *Hi/ 2 J =1/2 transitions were measured 
using 4 spectral windows each with 50 MHz of bandwidth 
sampled by 16384 channels. All of the 4.8 GHz measure- 
ments (0b — 153") utilized the position switching tech- 
nique with reference position located 30 arcmin west in 
azimuth from each source position. The H 2 CO 2n — 2i 2 
transition was measured using the dual-beam (9b = 51"; 
beam separation 330" in cross-elevation) receiver over 
50 MHz of bandwidth sampled by 16384 channels. The 
dual-beam system at 14.5 GHz allowed for both position 
switching and beam nodding. Both correlator configu- 
rations produced a spectral channel width of 3.052 kHz, 
which is approximately 0.2 and 0.08 km s _1 at 4.8 and 
14.5 GHz, respectively. 

To calibrate the intensity scale of our measurements, 
several corrections need to be considered: 

Opacity: Historical opacity estimates based on atmo- 
spheric model calculations using ambient pressure, 
temperature, and relative humidity measurements 
indicated that r at 4.8 and 14.5 GHz were ~ 0.01 
and 0.02 during our observations (assuming ele- 
vation > 30 degrees). The opacity corrections 
exp(r csc(EL)) are < 1.02 and < 1.04, respec- 
tively, and are considered negligible in our analysis. 

Flux: Assuming point-source emission, one can 
use the current relation (derived from 
measurements) for the aperture efficiency 

r\ A = 0.71exp(-(0.0163j/(Gffz)) 2 ) to con- 
vert antenna temperature to flux density. At 4.8 
and 14.5 GHz this yields r] A = 0.71 and 0.67, 
respectively. For elevation 90 degrees and zero 
atmospheric opacity, T^/S = 2.846?^ = 2.0 and 
1.9 for 4.8 and 14.5 GHz, respectively. These are 
the K/Jy calibration factors used to convert our 
spectra to flux assuming point-source emission. 
Note also that since the opacity correction is small, 
T'X = T A exp(Ar ) ~ 2k ~ T A , where m = 0.99 for 
the GBT. Using r\ m \, ~ 1. 327^, we can write the 
main beam brightness temperature as T m t, ~ 

The National Radio Astronomy Observatory is a facility of 
the National Science Foundation operated under cooperative agree- 
ment by Associated Universities, Inc. 
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TABLE 1 
Source List 



Source 


a 


<5 




T-i b 


T * dust 


Classification 




(J2000) 


(J2000) 


(km s ) 


(Mpc) 


( K ) 




NGC 253 


00-47-33 1 


—25:17:18 


251 


2.5 


34 


SAB(s)c 


NOC 520 


01 -24-35 3 


+03-47-37 


2281 


30.4 


38 




NOC fififl 


01 -43-01 7 


4-1 3-38-36 


856 


14.0 


37 


SRfsWnec HTT LINER 


NCC 891 


02-22-33 4 


4-42-20-57 


529 


10.0 


28 


SAfs'lh^-sn 

ja i o i u . ■ l> vj 


IC 342 


03-46-49 7 


-1-68-05-45 


31 


3.7 


30 


SARlVsVd 


Arp 55 


09-15-55 1 


-1-44-1 9-55 


11957 


159.4 


36 




wnr 9qd^ 

in vj\j ^yuo 






556 


6.2 


29 


SRl'tilH 

OJJ In iu 


UGC 5101 


09:35:51.6 


+61:21:11 


11810 


157.5 


36 


Syl.5, LINER 


M 82 


09:55:52.2 


+69:40:47 


203 


3.5 


15 


10, Sbrst 


NGC 3079 


10:01:57.8 


+55:40:47 


1150 


16.2 


32 


SB(s)c, LINER 


IR 10173+0828 


10:19:59.9 


+08:13:34 


14716 


196.2 




Sbrst 


NGC 3628 


11:20:17.2 


+13:35:20 


847 


9.6 


30 


Sb:pec:sp 


IC 860 


13:15:04.1 


+24:37:01 


3866 


51.5 




Sa, Sbrst 


M 83 


13:37:00.9 


-29:51:57 


518 


3.7 


31 


SAB(s)c 


IR 15107+0724 


15:13:13.1 


+07:13:27 


3897 


52.0 




Sbrst 


Arp 220 


15:34:57.1 


+23:30:11 


5434 


72.5 


44 


Pair, Sbrst 


NGC 6240 


16:52:59.0 


+02:24:02 


7339 


97.9 


41 


I0:pec, LINER, Sy2 


IR 17468+1320 


17:49:06.7 


+ 13:19:54 


4881 


65.1 




Sbrst 


NGC 6946 


20:34:52.3 


+60:09:14 


48 


5.5 


30 


SAB(rs)cd 



a Heliocentric velocity drawn from the literature. Source distance drawn from the literatu re (v < 
2000 km s — 1 ) or calculated assuming H = 75 km s^ 1 Mpc _1 . c From |Gao fc Solomon| j2004b| . d From 
NED a , Sbrst — starburst galaxy. 



"The NASA/IPAC Extragalactic Database (NED) is operated by the Jet Propulsion Laboratory, 
California Institute of Technology, under contract with the National Aeronautics and Space Ad- 
ministration. 



Source Structure: Several of our sources are known to 
have high-density structure (measured with high- 
dipole moment molecules like HCN) on scales ap- 
proaching the size of our beam at 14.5 GHz (51"). 
The beam coupling correction necessary to account 
for structure in our 14.5 GHz measurements, rela- 
tive to our point-source assumption, is given by: 



_ 01 + 01 

coupling ^ 



(i) 



This correction factor is less than 20% for 9s < 23". 
With the exception of M 82, none of the galax- 
ies in our sample have measured dense molecu- 
lar gas structure larger than ~ 20". We there- 
fore assume that, with the exception of M 82, all 
of the H2CO emission reported in these measure- 
ments is from structures smaller than the primary 
beams of our measurements. This assumption al- 
lows us to equate source main beam (T m {,) and 
radiation (Tr) temperature to our antenna tem- 
perature (T^) measurements. 

Absolute Amplitude Calibration: The GBT abso- 
lute amplitude calibration is reported to be 10-15% 
at all frequencies, limited mainly by temporal drifts 
in the noise diodes used as absolute amplitude cal- 
ibration standards. Relative calibration between 
our 4.8 and 14.5 GHz measurements is expected to 
be much better; estimated at ~ 5%. 

4. RESULTS 

In the following all spectra have been smoothed to 20 
km s" 1 (H 2 CO lio - In, OH, and Hllla) or 10 km s" 1 
(H2CO 2n — 212) to both increase the signal-to-noise ra- 
tio of our measurements and more closely match other 



mole cular spectral lin e measu rements of these galaxies 
(i.e. Gao & Solomon ( 2004a[ )). For each detection we 
list the peak intensity, velocity width (FWZI), and inte- 
grated intensity derived from direct channel-by-channel 
integration of the line profiles. 

4.1. H 2 CO 

Table|2]lists our H 2 CO lio — In and 2n — 2 i2 results. 
Spectra for the galaxies detected only in lio — In ab- 
sorption are displayed in Figure [2] while Figure [3] shows 
spectra from galaxies detected only in lio — In emission. 
IC 342, M 82 M 83, IR 15107+0724, and Arp 220 (Fig- 
ures [i] [5] [6[ [7] and U were all detected in both H 2 CO 
transitions. 

4.2. OH and Hllla 

Table [3] lists and Figure 9] shows our measured OH 
2 n 1/2 J = 1/2 F = 1 - andl - 1 hyperfine and Hllla 
radio recombination line results. For most of our mea- 
surements we did not detect the OH or Hllla transi- 
tions. The RMS for our OH 4765 and 4750 MHz and 
Hllla measurements of the undetected sources are listed 
in Table g] 

Rotationally excited OH has been used as a tracer of 



the m o lecular environmen t within AGN (Hcnkcl et al. 
jl986l , iHenkeh et ^jjj l987 ] ) |Henkel fc Wilson| ( [j^9rj| ~|rnT 



|pellizzeri et al. l 2006| )). AGJN come in two mam types 



those with (type 1J and those without (type 2) broad 
optical line emission. In the unified scheme of active 
galactic nuclei, all AGN are intrinsically similar. A sig- 
nificant column density of molecular material, in the form 
of a parsec-scale torus, obscures our view of the AGN in 
type 2 objects. Several surveys have attempted to detect 
the obscuring molecular material using measurements of 
molecular absorption or emissi o n, yielding few confirma - 
tions (e.g. Schmelz et al. ( 1986 1, Baan & Haschick ( 1992 1, 
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TABLE 2 
H2CO Measurements 



Source Transition T^/S peafc a v he ; a FWZI /T^dv//Sdv a 

(mK/mJy) (km s _1 ) (km s — 1 ) (mK km s — 1 /mjy km s -1 ) 



NGC 253 


lio 


- In 


-31.9(1.8)/-16.2(0.9) 


226.7(2.2) 


362 


-5782(165)/-2935(84) 


0.0074(4) 


NGC 520 


lio 
2n 


- In 

- 2i2 


-4.5(0.8)/-2. 3(0.4) 
(0.52)/(0.29) 


2287.6(6.5) 


207 


-469(61)/-238(31) 


0.0050(9) 
(0.0007) 


NGC 660 


lio 
2n 


- In 
-2i2 


-2.6(0.8)/-l. 3(0.4) 
(0.47)/(0.26) 


934.6(5.0) 


708 


-906(87)/-460(44) 


0.0025(8) 
(0.0006) 


NGC 891 


lio 


- Ill 


-2.7(0.8)/-1.4(0.4) 


499.5(5.0) 


323 


-445(71)/-226(36) 


0.0020(6) 
0.006(1)* 


IC 342 


lio 


- Ill 


-5.1(1.2)/-2.6(0.6) 


25.7(6.9) 


124 


-317(67)/-161(34) 




2u 


-2ia 


-2.9(0.6)/-l-5(0.3) 


26.7(2.8) 


104 


-140(16)/-78(9) 


0.0035(7) 


Arp 55 


lio 


- Ill 


(1.80)/(0.90) 








(0.002) 


NGC 2903 


lio 


- Ill 


(0.91)/(0.46) 








(0.0009) 


UGC 5101 


lio 


- Ill 


+3.5(0.8)/+l. 8(0.4) 


11864.6(6.3) 


99 


175(51)/89(26) 




M 82 


lio 


- Ill 


-23.0(1.0)/-11.7(0.5) 


214.3(5.0) 


430 


-4960(110)/-2518(56) 


0.0032(1) 
0.0015(2)* 




2n 


-2ia 


-5.2(0.7)/-2. 9(0.4) 


188.8(5.0) 


381 


-995(56)/-553(31) 


NGC 3079 


lio 


- Ill 


(1.30)/(0.66) 








(0.0009) 


IR 10173+0828 


lio 


- Ill 


(1.38)/(0.70) 








(0.002) 


NGC 3628 


lio 


- Ill 


-7.7(0.8)/-3. 9(0.4) 


840.8(5.0) 


368 


-1412(67)/-717(34) 


0.0070(7) 
0.0017(4)* 


M 83 


lio 


- Ill 


-2.4(0.6)/-l-2(0.3) 


522.9(10.8) 


155 


-183(42)/-93(21) 




2n 


-2i2 


-2.3(0.4)/-l-3(0.2) 


496.8(6.6) 


155 


-182(18)/-101(10) 


0.0026(5) 


IC 860 


lio 


- Ill 


+4.3(1.0)/+2.2(0.5) 


3892.5(8.1) 


267 


579(75)/294(38) 




IR 15107+0724 


lio 


- Ill 


+4.4(0.8)/+2. 2(0.4) 


3875.3(5.0) 


263 


569(73)/289(37) 






2n 


-2ia 


-3.2(0.5)/-l-8(0.3) 


3918.4(5.3) 


244 


-396(25)/-220(14) 


0.0041(6) 


Arp 220 


lio 


- Ill 


+5.1(0.6)/+2.6(0.3) 


5412.0(5.0) 


447 


1145(65)/581(33) 






2n 


-2ia 


-4.5(0.5)/-2.5(0.3) 


5379.2(4.6) 


410 


-923(34)/-513(19) 


0.0050(6) 


NGC 6240 


lio 


- Ill 


-2.6(0.4)/-1.3(0.2) 


7316.9(5.0) 


392 


-502(47)/-255(24) 


0.0026(4) 
(0.003)* 


IR 17468+1320 


lio 


- Ill 


(0.96)/(0.49) 








NGC 6946 


lio 
2n 


- Ill 

-2i2 


-2.4(0.6)/-1.2(0.3) 
(0.43)/(0.24) 


49.8(5.0) 


210 


-248(35)/-126(18) 


0.0022(6) 
(0.0005) 



a Heliocentric optical velocity frame, 
information. 



RMS noise levels arc for 20 (J = 1) and 10 (J = 2) km s 1 channels. Sec |5.2|for further 
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Fig. 2. — H2CO lio — 111 absorption spectra. The dotted line 
within each spectrum indicates the FWZI CO linewidth. 
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Fig. 3. — H2CO lio — 111 emission spectra. The dotted line 
within each spectrum indicates the FWZI CO linewidth. 



Staveley-Smith et al. (1992)). On the other hand, the 
existence of AGJN tori are con firmed by H2O megamaser 
emission in these objects (c/. Lo ( 2005| )). The observed 
OH line widths for the three sources detected in our sam- 
ple are similar, though slightly larger, than that mea- 
sured in the H2CO measurements of each source, sug- 
gesting a similar dynamical origin for the OH and H 2 CO 
emitting regions. Apparent optical depths for each of 
the OH and H2CO line region measurements range from 
0.04 (Arp 220) to 0.22 (IC 860), suggesting optically thin 
thermal absorption. 

4.3. Continuum Emission 

Table [5] lists our measured band-center continuum lev- 
els derived from the zero-level offset of our spectroscopic 
measurements. For most sources our measured contin- 
uum fluxes agree with those quoted in the NED archive. 
Exceptions to this agreement are IC 342: S4.8 = 0.087 Jy 
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TABLE 3 
OH and HI 11a Detections 



Source Transition T^/S peafe v kei FWZI /T^dv/JSdv 

(mK/mJy) (km s — - 1 ) (km s _1 ) (mK km s~ 1 /mjy km s _1 ) 



M 82 


Hllla 


29.0(1.4)/14.7(0.7) 


203.2(5.0) 


519 


7630(140)/3815(70) 


IC 860 


OH4750 


-3.6(0.8)/-1.8(0.4) 


3843.6(11.2) 


213 


-378(73)/-192(37) 




OH4765 


(0.77)/(0.39) 








IR 15107+0724 


OH4750 


-4.9(0. 6)/-2. 5(0.3) 


3880.2(13.8) 


410 


-1011(45)/-513(23) 




OH4765 


-2.2(0.6)/-l. 1(0.3) 


3891.6(13.0) 


351 


-380(45)/- 193(23) 


Arp 220 


OH4750 


-28.4(1. 2)/-14.4(0. 6) 


5407.3(3.1) 


486 


-6887(122)/-3496(62) 




OH4765 


-13.4(1.0)/-6.8(0.5) 


5449.9(4.5) 


575 


-3849(110)/-1954(56) 



Note. — Heliocentric optical velocity frame. RMS noise levels are for 20 km s 1 channels. 



TABLE 4 
OH and Hlllct Non-Detections 



Source 




RMS a ' b 






OH 4765 


OH 4750 


Hllla 




(mK/mJy) 


(mK/mJy) 


(mK/mJy) 


NGC 253 


2.22/1.11 


4.96/2.48 


5.56/2.78 


NGC 520 


0.76/0.38 


0.78/0.39 


0.78/0.39 


NGC 660 


0.88/0.44 


0.96/0.48 


0.78/0.39 


NGC 891 


0.92/0.46 


1.04/0.52 


0.84/0.42 


IC 342 


1.32/0.66 


1.28/0.64 


1.20/0.60 


Arp 55 


1.56/0.78 


1.62/0.81 


1.62/0.81 


NGC 2903 


1.12/0.56 


1.24/0.62 


1.24/0.62 


UGC 5101 


1.12/0.56 


1.08/0.54 


1.08/0.54 


M 82 


1.64/0.86 


1.46/0.73 


Emis 


NGC 3079 


1.46/0.73 


1.32/0.66 


1.34/0.62 


IR 10173+0828 


1.48/0.74 


1.82/0.91 


1.80/0.90 


NGC 3628 


1.06/0.53 


0.92/0.46 


0.82/0.41 


M 83 


0.78/0.39 


0.68/0.34 


0.86/0.43 


IC860 


0.78/0.39 


Abs 


0.96/0.48 


IR 15107+0724 


Abs 


Abs 


0.52/0.26 


Arp 220 


Abs 


Abs 


1.60/0.80 


NGC 6240 


0.72/0.36 


0.72/0.36 


0.84/0.42 


IR 17468+1320 


1.28/0.64 


1.60/0.80 


1.10/0.55 


NGC 6946 


0.60/0.30 


0.58/0.29 


0.74/0.37 


a RMS noise levels 


arc for 20 km s~ 


_1 channels. 


An "Abs" en- 



try indicates that the line was detected in absorption, "Emis 
in emission. 



(GBT) and 0.12,0.28 (NED), M 82: S i4 5 = 1.220/1.325 
Jy (GBT) and 1.79 Jy (140') or 1.7 Jy (NED), M 83: S 4 8 
= 0.347 Jy (GBT) and 0.7 Jy (NED), IC860: S 4 8 = 0.011 
Jy (GBT) and 0.040 Jy (NED), and IR17468+1320: S 4 8 
= 0.208 Jy (GBT) and 0.047 Jy (NED). For our opti- 
cal depth calculations, we use our measured GBT flux 
densities and assume an uncertainty of 15%. Continuum 
emission from starburst galaxies at these frequencies is 
likely due to synchrotron emission (S^ oc v~ a ) with spec- 
tral index a between 0.7 and 0.8. 

5. ANALYSIS 

5.1. Comparison to Previous Measurements 

NGC 253: Previously detected i n absorption in the 4.8 
and 14.5 GHz H2CO tr ansitions ( |Gardner fc Whiteoak 
1974| jBaan et ai||1990j ), but our GBT detection of the 
1 10 — In transition is signi ficantly different in shape from 
the line profile reported by Gardner & Whiteoak (19741. 
The interferometric mea surem ents of the lip — In tran- 
(19971 revealed an absorption 
size 39" x 12" in extent. 



sition by |Baan et al.\ 
structure with FWHM 
peak integrated flux measured by Baan et al. ( 1997 1 



The 



of 1.72 Jy beam -1 km s _1 is about 60% of the inte- 
grated flux we measure, suggesting that some extended 
structure detected in our GBT observations is missing 
in these interferometric measurements. Previous CS, 



H 2 CO, and high-J CO measurements (|Baan et al. |1990 
19971 iHiittemeister et a?.|| 1997| |Gusten et a/.||2006| esti- 



mate spatia l densities mthe r ange 10 4 b cm 



NGC 520: |Araya et aZ.| ( |2004| ) first detected the 1 10 — 1 41 
line in absorpti on, m good agreem ent with the GBT 
line properties 
NGC 520 at the 



Araya et al.\ (|2004| detected HI 10a in 
0.3 mJy level, but our GBT obser- 
vations of Hllla lack the sensitivity to detect a line 
of similar strength (0.25 mJy rms in 40 km s _1 chan- 
nels). Interferometric (6b — 6") mea surements of the CO 
J=l -> emission from NGC 520 ( |Sanders et~aD|1988 1 
revealed hig h concentrations of den se gas over a 12" (1.8 
kpc) region. |Yun fc Hibbard| ( |2001[ ) also imaged the CO 
J=l — > emission from JNGC 520 but with higher spatial 
resolution (Ob = 2". 4 x 2". 7), measuring dense gas over 
a somewhat smaller 7".l x 2". 9 region. 
NGC 660: Our GBT detection of the li - In line in 
absorption is marginal, but this line has previously been 
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TABLE 5 

Measured Continuum Levels at Band Centers 



Source 








14.5 GHz 






GBT 


Arecibo a 


NED 


GBT 


140' b 


NED 




(Jy) 


(Jy) 


(Jy) 


(Jv) 


(Jv) 


(Jy) 


NCIC 9^3 


1.795 




2.0 2.5 




0.37 




wr.r ^90 


083 


0.079 


0.1 


0.033 








0.153 


0.140 


0.19 


0.059 




< 0.09 


NOC! 891 


0.324 




252 342 








IC 342 


0.087 




1 2 28 


0.045 






Arp 55 


0.039 












wnr 9Q(i^ 
in uv ^yuo 


0.116 


0.041 


0.12 








UGC 5101 


0.061 




0.078 


0.032 






M 82 


3.266 




3.9 


1.220/1.325 


1.79 


1.7 


NGC 3079 


0.393 




0.32 


0.14 




IRAS 10173+0828 


-0.005 


0.020 










NGC 3628 


0.189 


0.131 


0.28 




0.11 




M 83 


0.326 




0.7 


0.075 






IC 860 


0.010 


0.016 


0.040 








IRAS 15107+0724 


0.019 


0.031 




0.022 






Arp 220 


0.213 


0.172 


0.2 


0.087 


0.11 




NGC 6240 


0.144 




0.178 








IRAS 17468+1320 


0.195 




0.047 








NGC 6946 


0.170 




0.2,0.5 


0.034 







The Areci bo 4.8 OH z continuum levels were obtained from spectral baselines by |Araya| 
i 2004) , Table [l"| b The NRAO 1 40' 14.5 GHz continuum levels were obtained from 



let al 



spectral baselines by |Baan et aZ.| l |1990[ l, Tabic [T] 




-1000 -500 500 

Velocity (km s ) 



1000 




Fig. 4. — H2CO211 — 2i2 (top) and lio — In (bottom) spectra of 
IC 342. The dotted line within each spectrum indicates the FWZI 
CO linewidth. 



-500 500 1000 

Velocity (km s ) 

Fig. 5. — H2CO 2n — 2i2 (top) and lio — In (bottom) spectra 
of M 82. The dotted line within each spectrum indicates the FWZI 
CO linewidth. 



et al. 


( 


1995 


iJ CO . 





CTTT^l 



Aalto 

0, J=2^ri7 

0; and HCN 



0, J=2 -» 1; C 18 J=l 
J=l — > emission from this galaxy and derived a mean 
source diameter of 14" from their CO J=l — > observa- 
tions. 



NGC 891: T his is a new detection of the lip — In line 
in absorption. Garcia-Burillo & Guelin (1995) measured 
the CO J=2 — I 1 and 1 — > emission toward this edge- 
on galaxy. The emission is extended along the galactic 
disk of extent > 120" x 20", but is mainly concentrated 
in a nuclear condensation ~ 20" x 10" in extent. 
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£ 




-500 500 1000 1500 

Velocity (km s ) 

Fig. 6. — H2CO 2n — 2i2 (top) and lio — In (bottom) spectra 
of M 83. The dotted line within each spectrum indicates the FWZI 
CO linewidth. 



£ 




4500 5000 5500 6000 

Velocity (km s ') 

Fig. 8. — H2CO 2n — 2i2 (top) and lio — In (bottom) spectra 
of Arp 220. The dotted line within each spectrum indicates the 
FWZI CO linewidth. 



£ 




3000 3500 4000 4500 

Velocity (km s ') 

Fig. 7. — H2CO 2n — 2x2 (top) an d lio — In (bottom) spectra 
of IR 15107+0724. The dotted line within each spectrum indicates 
the FWZI CO linewidth. 



IC 342: This is a new detection of H2CO: both lines 
are detected in absorption and show similar optical 
depth. Previous millimeter-wave length H2CO measure- 
ments (|Hiittemeister et al. 19971 estimate spati a l den 



cm 



Downes et al. 



(1992) 



sities in the range 10 4 
imaged the HCN J=l — > emission toward this galaxy, 
finding that the emission is composed of five condensa- 
tions with sizes ranging from 3.0" x 3.2" to 6.5" x 3.1". 
The total extent of these condensations is ~ 10" x 15". 
Arp 55: We did n ot detect H 2 CO , OH nor Hllla emis- 
sion in this galaxy. Sanders et al. ( 1988 ) imaged the CO 
J=l — > emission from this galaxy, detecting emission 
over a nuclear structure < 4" (< 2.8 kpc) in extent. 
NGC 2903: We did no t detect H 2 CO, OH , nor Hllla 
emission m this galaxy. |Helfer et al] (2003) imaged the 
CO 2 — > 1 emission from this galaxy! The CO emission 
emanates from a nuclear structure ~ 5" x 10" in extent, 
and from a more extended component ~ 3' x 10". 
UGC 5101: This is a new detection of the lio — In 
line in emission, and is the most distant H 2 CO emit- 
ter kn own, given the nondetections by Araya et al. 
( 2004} of previously claimed detections of H 2 CO in 



IKK5 12112+0305 and IRAS 10173+0828 (Baan et al. 



1993} 


UGC 


5101 


(IZhan 


g et al. 


2006 


HCN 


and HCO+ J 



Imanishi et al. ( |2006[ ) imaged the 
emission from this ULIRG, 
measuring a compact emission structure ~ 5" x 5" in 
extent. 
M 82: 

the lio - In l m e by Graha m et «Z. |( 1978| ) and in emission 

2 02 (218GHz) lines by |Baan| 



the 2 



H 2 CO was previously detected in absorption in 
line by Grah 
2i2 and 3 03 



11 



et al. (19901. We detect both the 4.8 and 14.5 GHz lines 



8 
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HI 1 la, OH4750, and 0H4765 




3000 3500 4000 4500 





IR 15107+ 


0724" 




0H4765 












3000 3500 4000 4500 



Velocity (km s 1 ) 

Fig. 9.— HIIIq, OH4750 MHz, and OH4765 MHz spectra. 



in absorption with high confidence (Figure [|). It is dif- 
ficult to concoct a physical scenario that would produce 
14.5 GHz emission at roughly tripl e the optica l depth 
of 4.8 GHz absorption as seen by Baan et al. ( 1990[ |, 
but our absorption line ratios are wholly compatible with 



LVG models (see ^5.3 1. The HllOa and Hllla lines, as 
well as other higher frequency hydrogen radio recombina- 



Bignell 1981 


Bell et al. 1984 


Rodriguez-Rico et a/.|2004 ) 


We re-detect the HI 11a line and confirm the two com- 



ponents identified by Bell et al. (1198 4 1. Previous CS and 



H^CO measurements (|Baan et at. 



990 



et «Z.||l997|) estimate spatial densities m t 



Hiittemeister 
re r ange 1 4 ~ " 



cm"". |Seaquist, Lee, fc Moriarty-Schieven (2006) im- 
aged the CO J =6 — > 5 emission from this galaxy, finding 
emission extended over a disk-like structure 40" x 15" in 
extent. 

NGC 3079: W e do not detect the li — In line detected 
in absorption by Baan et al. ( 1986 1 at —2.1 mJy; this line 
would appear as a '6.2a feature m the GBT spectrum, but 
no feature below —1.3 mJy (2. 0<r) is seen. Smoothing 
th e spectrum to ~ 10 km s _1 (the line width reported 



tion of the H 2CO lio — In line in emission, but |Araya| 



et al. ( 2004 1 do not detect the line at lower noise levels, 
and we confirm the non-d etection (but the GBT spec- 
trum would only detect the |Baan et al. ( 1993 ) line at the 
2.1cr level in ou r sm oothed 20 km s _i resolution spectra) . 
Planesas et al. (|1991[) imaged the CO J=l — * emission 



from this galaxy, detecting an unresolved nuclear source 
< 6" (< 6 kpc) in extent. 

NGC 3628: The lm — I n line was detecte d in absorp- 



tion by |Baan et al. ( 1986[ ) and confirmed by Araya et al. 
(|2004[ ), while the 2n — 2i2 line was not detected by Baan 



et al. 



In hue properties show 



fll990| ). The GBT 1 10 

good agreement with those measured by |Araya et al. 



fl2004| to within the calibration uncertainties. Stevens 
\et al.\ (|2005|) imaged the 450/um and 850/zm continuum 



emission from this galaxy and derived a core structure 
~ 30" in size. 

M 83: This is a new detection of H2CO: both lines are 
detected in absorption, and the 14.5 GHz line shows a 
significantly larger op tical depth than the 4.8 GHz line. 
|Sakamoto et al.\ ( |2004l ) imaged the CO J=2 -> 1 and 3 -> 
2 emission from this barred star burst galaxy, deriving a 
45" x 15" structure. 

IC 860: The lip — In emission line was detected 



by Baan et al. (19931, and the GBT detection shows 
good agreement in line properties (and somewhat worse 
agreem ent with the Arecibo observations of | Araya et al.\ 
( 2004 1) . We ar e not aware of any ima ging measurements 



Yao & Seaquist ( 2003[ ) measured the CO 
1 — > emission in 15 beams, which we 



of this galaxy. 
J=2 -> 1 and 

take as a maximum so urce size. 

IRAS 15107+0724: |Baan et al\ (11993b detected the 
lio — In line in emission, and the GBT line properties 
are in fair agreement with the discovery detection and in 
good agreement with Araya et a?.|(|2004l. We also detect 
the 2n — 2i2 line in absorption. |Planesas et al. ( 1991 1 
imaged the CO J=l — * emission from this galaxy, de- 
riving a compact nuclear source with size ~ 3". 
Arp 220/IC 4553: The lip - In line of H 2 CO was de- 
tected in emission by|Baan et a?.|fll98 6|, but the 2n — 2i 2 



line was not subsequently detected ( [Baan et al. 19901 



The sub-arcsecond resolution interferometric measure- 
ments of the lio — In transition by |Baan fe; Haschick 



( 1995 1 detected emission from two main components 
with individual sizes of ~ 0.5", ~ 1" in total extent, and 
integrated intensities of 83 and 682 mJy beam -1 km s _1 . 
The sum integrated intensity from these two components 
compares quite favorably to the total 1 10 — In integrated 
intensity and line width from our GBT measurements 
(581 ± 33 mJy km s _1 over a total line width of about 
500 km s -1 ). We also find good agreement between our 
4.8 and 14.5 GHz GBT spectra and those reported |Araya| 



et al.\ (120041). Recent CS, HCN, and HCO+ measure 



ments (|Greve et al. 20071 estimate spatial densities in 



the range 10 4 u cm 



The compendium of CO, HCN, HNC, HNCO, C S, and 
by [Baan et al. (1 1986 ) is 275 km s -1 ) yields a 0.4 mJy HCO+ mea surements of Arp 220 presented by |Greve 
LM S noise level 



RMS noise level and n o feature below . 5 mJ y at the 
expected line velocity. Stevens & Gear ( 2000 1 imaged 
the 450/im and 850/im emission from this galaxy, finding 
an unresolved core ~ 20" in extent. A similar source 



size was measured by Braine et al. (19971 in the 1.2mm 
continuum emission (p~ a 



17"). 



IRAS 10173+0828: Baan et al. (1993) claim a detec 



et al. ( 2007 1 point to an apparent dichotomy in the rel- 
ative strengths of the two velocity components exhib- 
ited by the spectra in these molecules. The two velocity 
components, at ~ 5350 km s -1 (western) and ~ 5500 
km s -1 (eastern) have been associated with tw o nuclei 
embedded within a circumnuclear r i ng or disk (|Scoville| 
et aZ.|1997||bownes fc Solomon|1998| . |Greve et «r] ( |2007| ) 
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point out that the western velocity component is stronger 
in the 12 CO J=2 -> 1 and HCN J=l -> transitions, 
but weaker in the HCN J=3 — > 2 and 4^3 transi- 
tions, than the eastern velocity component. In contrast, 
the eastern component is stronger in CS J=7 — ► 6 and 
HCO + J=4 — > 3 than the western component. Greve 



et al.\ (I2007J suggest that this dichotomy reflects differing 
spatial densities in these two components: the western 
component possessing lower spatial densities, the east- 
ern component possessing higher spatial densities. Our 
H2CO measurements do not necessarily reflect this trend. 
The lxo — In transition clearly shows a double-horned 
structure, but the 2 n — 2 12 transi tion does not (Figure [8]). 
NGC 6240: |Baan et al] ( |1993[) detected the 1 10 - In 
line in emission, but we detect the line in absorption. It 
is possible that spectral standing waves produced a spuri- 
ous emission line detection in the Baan et al. ( 1993) mea- 
surements. The GBT absorption line is detected with 
high confi dence (5.2a in a sin gle smoothed 20 km s~ 1 



channel). |Greve et al] ( |2007| noted that 12 CO, 13 CO, 



HCN, CS, and HCO" 1 " spectra all have a single gaussian 
profile which peaks at about the same velocity of 7160 
km s _1 . Our lxo - In spectrum peaks near 7300 km s _1 , 
slightly blueshiftcd relative to the syst emic velocity of 
7359 km s~ x . |Bryant fc Scoville| ( |1999| imaged the CO 
J=l — > emission from this galaxy and derived a source 
size of 10". 



IRAS 17468+1320: Baan et al.\ (|1993j claim a tenta- 
tive detection of the H2CO lxo — In line in absorption 
(3c). The GBT observations, with similar spectral res- 
olution, show no line, which would only appear at the 
2.5a level in our smoothed 20 km s" 1 spectra. We are 
not aware of any molecular spectral line or continuum 
emission imaging measurements of this galaxy. We arbi- 
trarily assume a source size of < 5" in our analysis. 
NGC 6946: This is a new dete ction of H 2 CO: th e 1 10 - 
lxi line appears in absorptio n 



Schinnerer et al. 
and iSchinnerer et al.\ (I2007I) imaged the CO J 
Land HCN 



(2006) 



J=2 — > 1, and HCJN J=l — > emission from this galaxy, 
detecting a compact nuclear source size of ~ 2" and a 
"nuclear spiral" structure 5" x 10" in size. 

5.2. H2CO Apparent Optical Depth Calculations 

The apparent peak optical depths (r, see Equation |2| 
for our H2CO absorption measurements (Table [2| have 
been calculated using the GBT continuum emission in- 
tensities listed in Tableland assuming T ex = 0. Cases 
where there are large discrepancies between the GBT 
continuum and other measurements are indicated by an 
asterisk in Table [2] We used the Arecibo continuum mea- 
surement for IRAS 10173+0828 (no continuum was de- 
tected by the GBT). Uncertainties in r reflect only the 
statistical error in the line measurement and neglect sys- 
tematic errors such as calibration uncertainties, but these 
should generally cancel out because the continuum and 
line measurements were obtained from the same spec- 
trum. The optical depth listed for H 2 CO non-detections 
in Table [2] is a la noise level. No optical depth is calcu- 
lated for emission lines. 

5.3. Spatial Density and Column Density Derivation 
Using LVG Models 

To derive the H 2 spatial density (number density) and 
H 2 CO column density of the dense gas in our galaxy sam- 



ple, we use a model which incorpora tes the Large Ve loc- 
ity Gradient (LVG) approximation (Sobolev ( 1960 1) to 
the radiative transfer. The detailed properties of our im- 
plementation of the LVG ap proximation are described in 
Mangum & Wootten (1993). This simplified solution to 
the radiative transfer equation allows for a calculation 
of the global dense gas properties in a range of e nviron- 
ments. As noted by Mangum & Wootten (1993), one of 
the major sources of uncertainty in an LVG model pre- 
diction of the physical conditions is the uncertai nty as- 
sociat ed with the collisional excitation rates used. IGreenl 
(1991) suggests that the total collisional excitation rate 
for a given H 2 CO transition is accurate to ~ 20%. There- 
fore, all physical conditions predicted by our LVG model 
are limited to an accuracy of no better than 20%. 

For the five galaxies where both the H 2 CO Lo — Li 
and 2xi — 2x2 transitions were detected we can derive a 
unique solution to the spatial density for an assumed gas 
kinetic temperature. By fitting the intercept between 
the H 2 CO Lo — In and 2xx — 2x2 transition ratio and 
the H2CO Lo — Li transition intensity (at an assumed 
gas kinetic temperature), one can derive a unique solu- 
tion to the (ensemble average) spatial density and H 2 CO 
column density, respectively, for each galaxy. For absorp- 
tion line measurements, the LVG model fitting procedure 
can also be done by using the measured apparent optical 
depths. For our galaxy sample, where T c < T cm f, (with 
one exception; M 82), this method has the disadvantage 
of involving an estimate of the excitation temperature 
(Tex ) in order to derive the apparent optical depth (see 
[ ]5.2| . Since T ex is dependent upon the input physical 
conditions, an accurate estimate of the apparent H2CO 
line optical depth is not possible without an assump- 
tion regarding the physical conditions in the gas being 
modelled. Therefore, we do not believe that fitting mea- 
sured apparent H2CO optical depth to those predicted 
by our LVG model is a better approach than fitting to 
measured H2CO line temperatures. In fact, as we show 
in £5.3.3 these two approaches to estimating the exci- 
tation of H2CO absorption lines in the presence of the 
cosmic background radiation and weak ambient contin- 
uum emission are roughly equivalent for the measure- 
ments presented in this paper. 

A model grid of predicted H 2 CO transition intensities 
over a range in spatial density, ortho-H2CO column den- 
sity per velocity gradient, and kinetic temperature of 



10 



4.5 



10 6 - 5 cm" 3 , 



10 



10.6 



10 



13.8 



cm 



7km 



and 



20-70 K, respectively, was calculated and compared to 
our H2CO Lo — Li and 2xi — 2x2 transition intensities 
(Table [2]). Table [6] lists the derived LVG model best-fit 
physical conditions for all of the galaxies in our sam- 
ple assuming T kinetic = T dufrf , with T dust from Table [T] 
T ex — 2.73 K, and negligible contribution du e to a ny 
background continuum emission (T c = 0). In { 5.3.3| we 
investigate how our results are affected by the addition 
of a background continuum source. 

A graphical representation of our LVG model deriva- 
tions of the spatial density and ortho-H 2 CO column den- 
sity is shown in Figure [TO] Note that: 

• The three galaxies with H2CO Lo — Li and 2xx — 
2i2 absorption (M 82, IC 342, and M 83) possess 
lower column densities and generally lower spatial 
densities. 
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TABLE 6 

Derived Physical Conditions 



Source log(N(ortho-H 2 CO)/Av) a log(n(H 2 )) T dust h 





(cm 2 /km s 1 ) 


(cm 3 ) 


(K) 


NGC 253 


11.98 + 0.02 




34 


NGC 520 


11.11 + 0.07 




38 


NGC 660 


10.88 + 0.10 




37 


NGC 891 


10.89 + 0.10 




28 


IC 342 


11.30+0.15 


5.05 ± 0.06 


30 


Arp 55 


< 11.20 




36 


NGC 2903 


< 10.90 




29 


UGC 5101 


12.30 + 0.10 




36 


M 82 


11.60 + 0.05 


4.69 + 0.03 


45 


M 82 (T c = 6.6 K) 


11.08 + 0.05 


4.90 + 0.03 


45 


NGC 3079 


< 11.06 




32 


IR 10173+0828 


< 11.08 






NGC 3628 


11.35 + 0.04 




30 


M 83 


11.45 + 0.24 


5.37 + 0.06 


31 


IC 860 


12.40+0.10 




10 c 


IR 15107+0724 


12.25 + 0.07 


5.68 + 0.01 




Arp 220 


12.48 + 0.01 


5.66 + 0.01 


11 


NGC 6240 


10.88 ± 0.10 




41 


IR 17468+1320 


< 10.95 






NGC 6946 


10.86 + 0.10 




30 



a For lio — In line only sources column density derived assuming n(H2) — 
10 J cm -3 (for absorption sources) or iq 5 65 cm -3 (for emission sources) and 
T K = 40 K. b From Tablc[l] c Assumed value. 



• The two galaxies with H 2 CO lio — In emission and 
2 U -2i2 absorption (IR 15107+0724 and Arp 220) 
possess the two largest spatial densities and column 
densities in our sample. 

Both of these features of our LVG model results are 
due to the unique spatial density dependence of the K- 
doublet transitions of H 2 CO. The higher density gas in 
IR 15107+0724 and Arp 220 is revealed by the H 2 CO 
lio ~ In emission. Furthermore, as pointed out in the 
individual source descriptions given in §5.1[ our derived 
spatial densities are well in-line with estimates made 
based on measurements done with other high-density 
tracers. 

As our spatial density and ortho-H 2 CO column den- 
sity calculations are based on an assumed kinetic tem- 
perature, an understanding as to how the spatial density 
and ortho-H 2 CO column density var y with changes in 
kinetic temperature is required (see S 5.3.1 L Following 



this analysis we investigate some nuances of our model 
fits for M 82 (potentially extended emission and strong 
background continuum) and Arp 220 (background con- 
tinuum driven lio — In maser emission) which require 
some additional discussion. 

5.3.1. LVG Model Dependence on Kinetic Temperature 

Using the same model fitting procedure described in 
§5.3| but allowing the kinetic temperature to vary from 
20 to 120 K, we can determine the uncertainty associated 
with our spatial density and ortho-H 2 CO column density 
calculations due to variations in kinetic temperature. If 
the actual kinetic temperature is 120 K rather than the 
assumed 40 K, for IC 342, M 82, M 83, IR 15107+0724, 
and Arp 220 the spatial density decreases by 30%, 30%, 
10%, 50%, and 40%, respectively. Similarly, a kinetic 
temperature of 120 K will cause the derived ortho-H 2 CO 
column densities to rise by 100%, 65%, 250%, 160%, and 



100%, respectively. To summarize, an increase by a fac- 
tor of three in kinetic temperature results in a 10-50% 
decrease in the derived spatial density and a correspond- 
ing 65-250% increase in the derived ortho-H 2 CO column 
density. 

5.3.2. Influence of Possible Extended Emission in M 82 



As revealed by CO 6^5 (Scaquist, Lee, & Moriarty- 


Schieven (2006b, 


Ward et al. ( 


2003}), HCiN 1 -» 


(Brouillet & Schil 


te (1993|)), anc 


I CS 1 -► and 2 -» 1 
32 possesses dense gas 


(Baan et al.\ 


1991 


]) images, M 



ing Equation [T] we calculate that the coupling correc- 
tion for this assumed emission distribution in our H 2 CO 
2n — 2i 2 measurements is /coupling — 1.2. This correc- 
tion would decrease the measured H 2 CO transition ratio 

( jTn(2ii-2i2)cfa ) k y 20% and produce a corresponding 
increase in the derived spatial density. Therefore, the 
derived spatial density for M 82 is a lower-limit to the 
true average spatial density in this galaxy. 

5.3.3. Influence of Background Continuum Emission 

The observed line temperature assuming perfect cou- 
pling of all continuum sources to the molecular cloud is 
given by 



T L = (T ex - T bg ) [1 - exp(-r)] 



(2) 



where Tt, g is the contribution due to all radiation sources, 
including the 2.7 K cosmic microwave background and 
any sources of continuum emission. By over an order- 
of-magnitude M 82 possesses the strongest continuum 
emission in our sample (Table pi), thus representing an 
extreme case in our analysis. We have modelled the in- 
fluence of a significant source of background continuum 
emission in our LVG models of M 82 by including an ad- 
ditional radiation source whose spectral characteristics 
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Fig. 10. — LVG model predictions for the spatial density (n(H2)) and ortho-H2CO column density per velocity gradient (N(ortho- 
H2CO)/Av) for the five galaxies in our sample with detected H2CO lio — In and 2n — 2±2- Two model fits are shown for M 82, 
representing models with no continuum and with T c = 6.6 K at 4.83 GHz (see text). Also shown are model line ratios (bold solid (positive) 
and dashed (negative) contours) and intensities (solid (positive) and dashed (negative) contours) for an assumed kinetic temperature of 40 
K and no background continuum emission (T c = 0). 



and intensities mimic our measured 4.8 and 14.5 GHz 
continuum emission. The effect that this additional ra- 
diation source has on our derived spatial density and 
H2CO column density is shown in Figure 10 Relative to 
the no-continuum model fit, the derived spatial density 
and H2CO column density are ~ 1.5 times larger and 
~ 3 times smaller, respectively, when the observed back- 
ground continuum emission is included. This is an ex- 
pected result, as the relatively stronger continuum emis- 
sion source at 4.8 GHz provides a larger contribution 
to the excitation of the observed H2CO lio — In line 
emission, lowering the model predicted transition ratio 

( jTH(2i°-2i2)<fo )" Therefore, ignoring contributions due 
to background continuum emission produces a derived 
spatial density for M 82 which is a lower-limit to the 
true average spatial density in this galaxy. 

A further check of the influence of the continuum emis- 
sion from our galaxy sample in the excitation of the 



H2CO molecule is afforded by a comparison between 
the apparent H 2 CO optical depth estimates and those 
predicted by our LVG models where we have assumed 
T c = 0. For the three galaxies where both the lio — In 
and 2 11 — 2 12 transitions are observed to be in absorp- 
tion (IC 342, M 82 with T c = 6.6 K, M 83), line optical 
depths, within both the modelling and optical depth cal- 
culation uncertainties, are consistent. For the two galax- 
ies where only the 2n — 2 12 transition is observed in ab- 
sorption (IR 15107 and Arp 220), LVG model-predicted 
optical depths are a factor of two larger than the calcu- 
lated apparent optical depths. Since the apparent optical 
depth calculation assumes that T ex — 0, our apparent op- 
tical depth estimates are likely to be lower limits to the 
true line optical depth. Given the approximations made 
to derive the apparent and LVG modelled optical depths, 
we suggest that there is good consistency between our 
continuum-less model-predicted optical depths and those 
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estimated by our apparent absorption line optical depth 
calculations. 

5.3.4. H2CO Excitation Pumping by Continuum Emission 

in Arp 220 

Bright and compact nonthermal continuum emission 
at frequencies near 4.8 GHz can in duce maser excitatio n 
of the H 2 CO 1 10 - In transition ( [Baan et al] Ql986| )). 
This maser excitation would amplify the background 
continuum emission and be observed as anomalously 
bright H 2 CO (maser) emission. Assuming an H 2 CO and 
nonthermal continuum emission source size of ~ 0.25" 
and nonthe rmal continuum brightness temperature of 
- 2xl0 5 K, |Baan et JT] ( |1986[ ) derived an H 2 CO li -lu 
brightness temperature of ~ 3 x 10 4 K. To study the 
conditions under which such a high brightness tempera- 
ture would be produced Baan et al. constructed an LVG 
model which included the effects of a strong nonthermal 
radiation field. These models showed that weak ampli- 
fication (t4.sghz > —0.025) of the lio — In transition 
results for n(H 2 ) < 10 5 cm -3 and H 2 CO abundances of 
> 10- 9 . 

Our measurement of H 2 CO 2n — 2i 2 absorption in 
Arp 220 puts some significant constraints on the poten- 
tial for maser amplification of the lm — In transition. We 
ha ve reproduced the LVG model predictions described 



which lead to amplification of the lio — In transition 
also predict: 

• Weak (ti4 5 qhz — 0.5 x T4.sgHz) maser amplifica- 
tion of the 2n — 2i2 transition, 

transition intensity ratios which are not 



2n — 2i 

consistent with our measured results. 



This result forces us to conclude that either: (1) the 
lio — In transition is not masing or (2) the two K-doublet 
transitions do not trace the same dense gas component. 
Note that this result is i ndepe ndent of the assumed ki- 
netic temperature (see S; 5.3.1 1. Without high spatial 
resolution measurements oi the 2n — 2i 2 transition we 
cannot discount the non-cospatiality conclusion. Based 
on the nearly identical line profiles from 5200 to 5700 
km s -1 , we believe that the two K-doublet transitions 
trace the same dense gas component. Therefore, we con- 
clude that the lio — In emission in Arp 220 is not due 
to maser emission, but does represent the signature of 
a dense gas component near the two starburst nuclei in 
this galaxy. 



5.3.5. Single Transition H2CO Column Density Estimates 

For a given spatial density and kinetic temperature a 
single H 2 CO transition can be used to derive the H 2 CO 
column density. To estimate H 2 CO column densities 
toward the galaxies that we have only detected in the 
lio — In transition, we assume a spatial density of 10 5 
cm~ 3 (for absorption sources) or 10 5 - 65 cm -3 (for emis- 
sion sources) and kinetic temperature of 40 K. Our spa- 
tial density assumption is based upon our derivations 
of the spatial density in the five galaxies where both 
H 2 CO transitions were measured, while our assumed ki- 
netic temperature is typical of our galaxy sample (see 
Table [T). Using our measured lio — In intensities we 



- log(n(H 2 )) 


= 5.0 cm 3 


^—i — * 


T K = 40 K 




NGC 253 




— 1 NGC 3628 




NGC 520 t""F 


NGC 891 

NGC 660 
NGC 6240 
NGC 69-16 





-0.02 

< (1,0- In) (K) 



Fig. 11. — LVG model predictions for the ortho-H^CO column 
density per velocity gradient (N(ortho-H2CO)/Av) for the seven 
(large font) galaxies in our sample with detected H2CO lio — 111 
absorption only. We also show upper limits (3<r) to N(ortho- 
H2CO)/Av for our five non-detected (small font) galaxies. Not 
shown on this plot are the two lio ~ 111 emission line sources 
IC 860 and UGC 05101. 



Baan et al. (1986), but also find that the conditions ure 



derive the H 2 CO column densities for each galaxy (Fig- 
and Table [§. 



11 



5.4. A Cautionary Note Regarding Critical Density 

Very often estimates of the spatial density in molecular 
clouds are based on a calculation of the "critical density" 
of a measured molecular transition. The critical density 
of H 2 (n crit ) is given by: 



n crit (cm ) 



Auiis- 1 ) 



J2i C ui (cm 3 s x ) 



(3) 



where A u i is the spontaneous emission (Einstein A) rate 
for the transition between levels u and 1, while C U i is 
the H 2 collisional excitation rate out of the upper energy 
level into level i. This critical density is really a measure 
of the relative importance of radiative optically-thin and 
collisional excitation to the population of the upper state 
in a molecular transition. A u \ is given by: 



1{S )= o |Wt 



3/lC 3 

64tt 4 i/ 3 
3/ic 3 



S11 2 



(4) 
(5) 



where |/Lti u | 2 is the dipole moment matrix element for the 
transition invo lving levels 1 and u, [x i s the dipole moment 
(2.331 Debye; |Kondo fc Oka] ( [i960] )), and S is the line 
strength, given by: 



S 



K 



J(J + 1) 



for (J,K) -> (J,K) 



(G) 



where we use the (very good) approximation that H 2 CO 
is a symmetric top molecule. For the lio — 111 & n d 
2n — 2i2 transitions S = | and |, respectively. Using 
these line strengths, along with [i, in Equation [5] yields 
A ul = 3.6 x 10~ 9 and 3.2 x 10" 8 s~\ respectively. The 
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sum over the collisional excitation rates 2 out of level u 
(EiC«i) are 4.4 x 1CT 10 and 4.2 x 1CT 10 cm 3 s"\ re- 
spectively, at a kinetic temperature of 40 K. This implies 
that the critical densities for the lio — In and 2n — 2%2 
transitions are 8 and 76 cm -3 , respectively. These criti- 
cal densities arc at least three orders-of-magnitude lower 
than critical densiti es for low-lyin g: transi t ions o f HCN. 
HCO+, and CS (c/. ~ 



Janscn 



dl995|,|Evansl(|l999|) which 



have been used as high-density tracers of star formation 
As we are clearly measuring densities which are many 
orders of magnitude larger than the critical density for 
the H2CO transitions studied in this work, the critical 
density appears to be a poor indicator of the density 
probing properties of a molecular transition. Since sig- 
nificant optical depth in a transition can lower the critical 
density by limiting the radiative contribution to the up- 
per energy level population, the critical density is not 
actually a lower limit to the density required to excite a 
molecular transition under many interstellar conditions. 
Clearly density estimates based on critical density argu- 
ments are nothing more than crude estimates to the true 
spatial density in a molecular cloud. 

5.5. H2CO Luminosity and Dense Gas Mass 

Studies of the correlation between the infrared and 
molecular line luminosity in luminous infrared galax- 
ies have been used to characterize the power source for 
their intense infrared emission. Surveys of the HCN 



( Gao fc Solomon! (l2004b|),|Gao fc So lomon (|2004bJ) and 



HCO 



(|Garcia-Carpio et al. 
these galaxies have derived a tight correlation between 
their molecular line and infrared luminosities. A com- 
parison between the infrared and HCN J=l — > emis- 
sion in our Galaxy has extended this correlation down 
to Galactic dense core scales, extend ing the correlatio n 



-C 



(2006)) 



emission toward 



over 7-8 orders of magnitude in Lr/j ( |Wu et al. | (2005 1) 
This correlation has been interpreted as evidence that 
star formation is the main power source for the large 
infrared luminosities observed in these galaxies. 

Following the formalism for calculating molecular 
spectr al line luminosity elucidated by |Solomon et al.\ 
(19971, we note that the monochromatic luminosity 
,j(is res t), the observed flux density S(u obs ), and the lu- 
minosity distance Dl are related by v rest L(u rest ) = 
^ S(v obs )v obs D 2 L . Relating L{v rest ) to the spectral line 
luminosity L Une : 



)dv 



-A-KSu ne (v bs)VobsD 2 L 



dv 



1.043 x lO- 3 S lme Avv rest (l + z)~ x D\ L G (7) 



with Sune in Jy, Av in k m s 
and Dj} is in Mpc. As 



'rest 
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V nhs(l+ z) in GHz, 



(1997) point out, 



it is convenient for relating spectral line integrated inten- 
sity measurements to introduce the quantity L' , which is 

2 Note that H2CO collisional excitation rates are calculated us- 
ing He as collider. The quoted rates have been scaled by a factor 
of 2.2 to cor rect for the difference between He and H2 collisions 
( |Green|199l" i. 



J = Da (1 + z ) 2 i where D^ is the luminosity distance and 
Da is the angular size distance. 



the source-integrated surface brightness. Two spectral 
lines with the same main beam brightness temperature 
(T mb ) and spatial extent will have the same L' , making 
the L' ratio of two spectral line measurements an indica- 
tor of the physical conditions in the gas. 

For any molecular spectral line measurement L' is 
given by: 

L' = I T h dvdtt x D 2 



T mb dvQ, s * b (l + z)D 2 A 
23.504n„ b (arcsec 2 )(l + z)- 3 D 2 L (Mpc) 
T mb (line)dv K km s~ l pc 2 



_ 23.504 tt(6> 2 ( arcsec) + 9 2 B (arcsec)) 3 
= 41n(2) + ° 

D 2 L (Mpc) J T mb (line)dv K km s" 1 pc 2 (8) 
= + / S lm e dv 

lKV rest J 

= 3.256 x lOV-yi + z)- 1 ^ 

with S Hne in Jy, T mb (line) in K, dv in km s 1 , v rest = 
v bs{^+z) in GHz, Dl in Mpc, and fl s * b the solid angle of 
a gaussian source convolved with the gaussian telescope 
beam, given by: 



9. 



s*b 



ne% (e 2 s + e 2 B 



41n(2) 



I 2 



(10) 



Studies which compare with L' are seeking to in- 
vestigate the correlation between source-integrated dense 
gas column density and mass to IR luminosity. Since we 
have LVG-derived H2CO column densities, the mass of 
dense gas traced by H 2 CO can be derived as follows: 



m -no n 2 ^ TOg 2 

lvl dense — I ^mol il S' U A v 



: 23.504fi s (arcsec 2 )(l + z)- 2 D L (Mpc) 

N mo i(cm~ 2 )ixm H2 {gm) 2 _ 2 

pc cm 

1.125 x 10 5 tt9 2 (arcsec) D 2 (Mpc) 



4ln(2)(l + z) 2 
Nmoi (cm~ 2 ) fj,m H2 (gm) 



Mr. 



(11) 



where is the mass of molecular hydrogen, p accounts 
for the mass fraction due to He, and X mo i is the abun- 
dance (relative to H 2 ) of the molecule. 
For emulateapj 

In Table [7] we list the infrared luminosities derived 
from measurements culled from the literature along with 
the dense gas mass (M ( iense) derived from HCN and our 
H 2 CO measurements. The infrared luminosity Ljji over 
the wavelength range from 8 t o 1000 /xm is derived from 
the Sanders & Mirabel ( 1996 1 relation 
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TABLE 7 

Infrared Luminosities and Dense Gas Masses 



Source 


Size a 


T b 

Lib 


M dense (HCN) c M dense (H 2 CO) 




(arcsec) 


(1O 1O L ) 


(10 8 M(7i) 


(10 8 M o ) 


NGC 253 


39 X 12 


2.78+0.001 


2.7 


3.02+2.14 


NGC 520 


12 


8.10+0.05 




10.16+7.37 


NGC 660 


14 


3.10+0.01 


> 2.6 


5.84+4.34 


NGC 891 


20 x 10 


1.86+0.005 


2.5 


1.39+1.04 


IC 342 


10 x 15 


1.47+0.002 


1.7 


0.14+0.11 


Arp 55 


< 4 


45.37+1.33 


38 


<56.91 


NGC 2903 


~ 5 x 15 


1.55+0.004 


> 0.9 


<0.22 


UGC 5101 


~ 5 


83.11+1.03 


100 


16.61+12.25 


M 82 


40 x 15 


5.92+0.001 


3.0 


3.81+2.73 


NGC 3079 


~ 20 


5.37+0.02 


~ 10 


<25.27 


IRAS 10173+0828 


< 6 


57.17+1.79 




<485.44 d 


NGC 3628 


~ 30 


1.78+0.008 


2.4 


19.56+13.95 


M 83 


45 x 15 


1.26+0.001 


3.5 


1.14+1.05 


IC 860 


< 15 


14.48+0.10 




< 581.19 


IRAS15107+0724 


~ 3 


18.25+0.10 




32.22+23.37 


Arp 220 


~ 1 


155.76+0.32 


30-370 


19.85+14.05 


NGC 6240 


10 


67.99+0.46 


80-260 


77.39+57.60 


IRAS 17468+1320 


< 5 


8.71+0.18 




<5.93 e 


NGC 6946 


5 x 10 


1.43+0.002 


4.9 


0.06+0.05 



a See j j5.l| for re ferences. 13 IRAS flux es fro m|Sanders et al.\ho03\. c Except for Arp 220 
and NGC 624 l|Greve et al. [2007 , from |Gao fe Solomon| j2004a| l d Assumed H 2 CO 
lincwidth from |Araya et aL| l |2004p Assumed H 2 CO linewidth from |Baan et "oI1 l |1986[ l 



L IR = 4TrD 2 L F IR (12) 
F IR = 1.8x 10- 14 (13.48/ 12 

+5.16/ 2 5 + 2.58/ 6 o + hm)W m - 2 (13) 

with /12, / 2 5, jkb /ioo given by the IRAS 12 /im, 25 /xm, 
60 /im, and 100 ^m fluxes in Jy, respectively. The H 2 CO- 
derived dense gas masses were calculated using Equa- 
tion ITTI with N(ortho-H 2 CO)/Av from our LVG model 
fits (Table k3| , Av derived from our direct-measurements 
(Table [§, and assuming X(H 2 CO) = 10" 9 . In addi- 
tion to the derived uncertainties in N(ortho-H 2 CO)/Av 
and Av, we have conservatively applied 50% uncertain- 
ties each to our source sizes and H 2 CO abundances when 
calculating M dense- Due to the lack of H 2 CO source 
structure information our H 2 CO-derived Mdense values 
are only estimates to be compared with other molecule- 
derived dense gas masses. 

Our Mdense estimates range from 6 x 1O 6 M 
(NGC 6946) to 8 x 1O 9 M (Arp 220). In general, M dense 
derived from our H 2 CO measurements agrees to within a 
factor of 3 with Mdense derived from HCN measurements. 
In a few cases, such as IC 342, UGC 05101, NGC 3628, 
and NGC 6946, the dense gas mass derived from H 2 CO 
is as much as a factor of 30 too low relative to Mdense 
derived from HCN. In the case of IC 342 we can also 
calculate Mdense using our derived n(H 2 ) (Table [6j. As- 
suming a thin-disk for the emitting volume the mass is 
derived from 

^fl2 



ML 



-Dfhm H2 n(H 2 ) 



(14) 



41n(2)' 

where h is the thickness of the galactic disk. For IC 342 
we find that M' dense ~ 2 x 10 9 Af© assuming ft, = 20 pc 
(~ \8s)- This value for the dense gas mass is more in- 
line with other estimates of Mdense (i-e. from HCN). One 
way to bring our estimates of Mdense and M' dense into 
agreement is to use a lower H 2 CO abundance of 10~ 10 . 



Gao fc Solomon| ( |2004a) found that M de „ se (HCN) are 
factors of 5-200 smaller than M(iJ 2 ) derived from CO 
measurements for their sample of spiral, LIRG, and 
ULIRG galaxies, many o f which are included in ou r sam 
pie. As pointed ou t by Gao fc Solomon (2004a), Gac 
fc Solomon] ( |2004b[ ), and |Greve et al.\ ( |2007| ), the pro- 
gressive difference between the mass traced by dense gas 
tracers such as CO, HCN, and H 2 CO reflects the hier- 
archical structure of the giant molecular clouds in these 
galaxies. Our Mdense (H 2 CO) values suggest that H 2 CO, 
like HCN, traces a denser, more compact, component 
of the giant molecular clouds in our galaxy sample than 
CO. This result is consistent with high-resolution studies 
of the K-doublet H 2 C O emission in our own Galaxy (cf. 



Mangum et al. 



(1993)). 

show 



In Figure 12 we show Mdense ver sus Lj R . Note 
the tre nd, similar to that derived by |Gao fc Solomon] 



(2004a) using HCN J=l — > measurements. We also 
show in Figure 12 the log-linear relation log(L/#) = 
2 + l.Qlog(M ( j ense ). Note that given the uncertainties 
regarding dense gas mass source size we have not per- 
formed a formal fit of this relation. The linear correlation 
between the dense gas mass traced by H 2 CO and the in- 
frared lu minosity is consistent with the c orrela tion first 
noted by Solomon, Downes, fc Radford| ( 1992 1. Active 
star formation in 1R- bright galaxies is clearly driven by 
the amount of material available to form stars. Further- 
more, the relationship to Kcnnicutt-Schmidt laws, which 
relate the star format ion rate t o the mass of gas avail- 
able to produce stars (Schmidt (1959 1963[ )) in galaxies, 
to these measurements becomes apparent in the Lj R -to- 
M d 



son 



correlation. As shown by Krumho lz fc Thomp-| 



(12007]), the actual form for this correlation depends 
upon the spectral line chosen to derive the spectral line 
luminosity L mol . Since M dense cx L moh the suggestion 
of a linear correlation between L\ R and MdenseifliCC?) 
implies that, similar to HCN, H 2 CO traces the dense 
star-forming gas in starburst galaxies. 
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Fig. 12. — Ljr versus M^ ense . Mj ense has been derived using Equation 1 1 1 1 assuming X(H2CO) = 10~ <J and uncertainties calculated 
assuming 50% errors in Cl a and X(H'2CO). Filled circles and squares are used to indicate H2CO-detected galaxies at distances > and < 
50 Mpc, respectively. The dotted line represents the log-linear relation log(Ljn) = 2 + 1.0 log(M ( ; ense ). 



A related trend can be seen when one compares Ljr to 
our derived n(H2) for the five galaxies within which we 
have derived spatial densities. From Tables [6] and [7] we 
note that the three galaxies with Ljr ~ 1 — 5 x 10 L@ 
have derived spatial densities (n(i?2) — 10 5 cm~ 3 ). The 
two higher luminosity galaxies Lir > lO 11 !/© have 
higher derived spatial densities of n(H 2 ) = 10 5 7 cm -3 . 
Even though our number statistics are small (five galax- 
ies), there appears to be a trend toward higher spatial 
density for galaxies with higher infrared luminosity. This 
is likely another representation of the Lj^-M^ense corre- 
lation. 

6. CONCLUSIONS 

Using measurements of the 1 10 — In and 2n — 2 12 K- 
doublet transitions of H2CO we have derived accurate 
measurements of the spatial density (n(H2)) in a sample 
of starburst galaxies. The derived densities range from 
10 4 7 -10 5 7 cm~ 3 , consistent with the suggestion that the 
high infrared brightness of these galaxies is driven by ex- 
treme star formation activity. We believe that these spa- 
tial density measurements are the most accurate mea- 
surements of this important physical quantity in star- 
burst galaxies made to-date. We have also used our 



H2CO measurements to derive a measure of the dense 
gas mass which ranges from 0.06 — 77 x 10 s Mq, gener- 
ally consistent with previous measurements, mainly from 
studies of the HCN emission in these galaxies. The lin- 
ear correlation between the IR luminosity (Lm) and the 
dense gas mass (Mdense) found in larger samples of the 
HCN emission in starburst galaxies is also apparent in 
our H2CO measurements. This further supports the sug- 
gestion that active star formation in IR-bright galaxies 
is driven by the amount of material available to form 
stars. We also note a related trend between and 
our derived n(i? 2 ) for the five galaxies within which we 
have derived spatial densities. The three galaxies with 
lower IR luminosities have lower (n(i? 2 ) — 10 5 cm~ 3 ) 
derived spatial densities, while the two higher luminos- 
ity galaxies have higher (n^H^) — 10 5 7 cm" 3 ) derived 
spatial densities. This is likely another representation of 
the LjR-Mdense correlation. 
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